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The initial production and dynamical expansion of hot spherical nuclei are ex- 
amined as the first stage both in the projectile-multifragmcntation and in central 
collision processes. The initial temperatures, which are necessary for entering the 
adiabatic spinodal region, as well as the minimum temperatures and densities, 
which are reached in the expansion, significantly differ for hard and soft equations 
of state. Additional initial compression, occurring in central collisions leads most 
likely to a qualitatively different multifragmentation mechanism. Recent experi- 
mental data are discussed in relation to the results of the proposed model. 



1 Introduction 

Spinodal instabilities as a possible mechanism of multifragmentation was first 
suggested by Bertsch, Skinens and CugnonEla. Recent data from the ALADIN 
and FOPI experiments ErQuQ shed some light an such a possible explanation. 
In particular the discovery of the caloric curveB, which looks like the one for 
a liquid-gas phase transition in water, suggests to study the mechanism of 
spinodal decomposition in more detail. 

The expansion of the excited (projectile) residue can be regarded as the 
first stage in the multifragmentation process. For large enough initial excita- 
tion energies and temperatures the nuclear matter of the residue reaches states 
of mechanical volume instability, where small density fluctuations grow expo- 
nentialliMtj. For an illustrative analysis of the instability evolution we refer 
to ref. Ej. J^he .fragmentation process itselLhas been studied within molecular 
dynamics EiJ£3 : E3Ej. Studies by Friedman E£rt3 within an expanding emitting 
source model (EES) show that intermediate mass fragments (IMFs) are in- 
deed created within time intervals of about 50 fm/c indicating a simultaneous 
breakup of the residue. 

We have been investigating, already for some time now E3, the expan- 
sion dynamics of spherical hot nuclei for peripheral collisions. For low initial 
excitation energies one obtains monopole vibration with the largest possible 
amplitudes. Around an initial temperature of 5 MeV (i.e. E* /A = 2 MeV) the 
matter inside the expanding drop reaches, e.g. in the T, g plane, the bound- 
ary of the spinodal regime. The important feature at higher initial excitation 
energies is the occurrence of turning points deep in the spinodal region. The 
system stays long enough around these turning points, which all are located 
around T = 4 MeV, such that instabilities can develop and fragments can be 
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Figure 1: Time evolution of the Au residue for a soft equation of state in a peripheral 
collision (without initial compression) in the temperature — density plane. The crosses 
on the trajectories indicate time steps of 5 fm/c. The shadowed zone corresponds to the 
adiabatic mechanical instability region for infinite nuclear matter. The numbers give the 
charge numbers initially and at the turning points which am connected by the heavy solid 
line. The open boxes denote the experimental data El (without error bars). 

formed. At initial temperatures T > 12 MeV the system runs through the 
freeze-out density of 0.3go before reaching a turning point. 



2 The Model 

In this section we describe a few essential features of the pjodel. A more 
detailed description of the model will be published elsewhere i3. 

First we study the peripheral collision setup. We start out from a ther- 
malizcd spherically symmetric projectile-like nucleus at ground state density 
which has been creAtedJii the peripheral high-energy collision according to 
the abrasion modelEjcErEll. The typical temperature for such a hot nucleus is 
not too high, such that the mean free_Qath of nucleons is still comparable or 
even larger than the size of the systemtJ. Thus, we can approximately assume 
homogeneity in density and temperature for the expanding nucleus. 

During the expansion the hot nucleus evaporates particles. We consider 
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only neutron and proton evaporation. Dcutcrons may be considered to arise 
from coalescence of nucleons. We follow the evolution of the excited remnant 
and plot the trajectories in the temperature - density plane. 

If the excitation energy is high enough, the expansion leads into the in- 
stability region of the equation of state. For not too high excitation energies 
the expansion of the homogeneous system stops at a certain density. Around 
these turning points the system has enough time (> 30 fm/c, Era) to develop 
inhomogeneities in density, and subsequently decay into many fragments. The 
final mass and charge spectra is not given by this model, however some quan- 
titative comparison with experiments concerning temperature and excitation 
energy (caloric curve) can be made. 

For a central collision setup we considered a somewhat similar scenario. 
Here the initial expanding source can achieve much higher initial temperatures, 
and also one should take into account conxoression. The initial state was 
estimated on the base of the shock model E3. Due to the compression the 
evolution of such a system is so fast, that the evaporation can be neglected 
(the "evaporated" particles cannot leave the excited source), and except very 
low excitation energies no turning is obtained. The calculation is stop at the 
break-up density of 0.3po obtained from analyzing the ALADIN data. 

3 Results 

According to the abrasion modeB, the initial excitation energy (temper- 
ature) is correlated with the charge number of the projectile residue. These 
numbers are given on the right and upper parts of Fig. |l|. The calculated ex- 
pansion trajectories in the (T, g)-plane are significantly different for a soft and 
a stiff equation of state (EOS). 

We follow the trajectories up to their turning points if reached within 200 
fm/c. The turning points are indicated by the heavy solid lines. Turning 
points exist for initial temperatures up to about 14 MeV. We notice that the 
turning points for the same initial excitation lie one to two MeV in temperature 
and O.Iqo to 0.2qo in density higher for the stiff EOS as compared to the soft 
one. Moreover, we see that the turning points are located almost at the same 
temperatures independent of the initial excitation, i.e. around 4 MeV and 5 
MeV for the soft and hard EOS, respectively. This plateau does not move if 
one changes the initial excitation energy, the only effect is, that the charge 
numbers become different initially and at the turning points. 

A natural criterion for multifragmentation to occur, is whether the system 
reaches the region of volume instability (adiabatic spinodal region) , where the 
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Figure 2: The caloric curve as deduced by the ALADIN collaboration^ (crosses) compared 
with the pre-cluster condition. The total excitation energy includes flow energy which is 
significant in the isentropic dynamics for E* /A > 15 MeV. The discrepancy at large excita- 
tion energies may therefore indicate that dissipation is not negligible during the expansion 

process. 



derivative of the pressure with respect to the density along an adiabat becomes 
negative. Since the system is closed we consider the adiabatic process to be 
the relevant one. In Fig. the region of instability is shadowed and bounded 
by the spinodal. To enter the spinodal region and stay there for more than 30 
fm/c one needs initial temperatures of about 8 MeV and 11 MeV for the soft 
and hard EOS, respectively. 

Because of the occurrence of turning points in the expansion, the mul- 
tifragmentation process of the projectile residue has a unique feature: the 
subsequent decay into fragments is expected to be rather free from collective 
flow. This is different from the fragmentation of compressed compound sys- 
tems which are formed in central collisions. As our calculations show already 
for a small initial compression (1.5 go) no turning points occur for realistic ini- 
tial excitation energies. Thus the fragmentation process becomes qualitatively 
different and looks more like an explosions. 

In Fig. H we compare our results for the temperature and the total excita- 
tion energy in the remnant with experimental data from the ALADIN exper- 
iment. Our points are determined at the turning points if the corresponding 
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Figure 3: The narloric curve as deduced from the FOPI experiments at 100, 150 and 250 
MeV/u collisionsa'EI (crosses) compared with our results (circles for Qt, r =0.3 go and diamonds 
for g(, r =0.25 go)- The flow energy is subtracted from the excitation energy. 



density is higher than Qb r ~ 0.3go> or at the break-up density otherwise. For 
very low excitation energies the average values of the vibration are chosen. 
There are some deviation at high excitation energies, which is due to the aris- 
ing flow energy in our model. 

In Fig. |^ the comparison with the FOPI data uB are plotted. The ex- 
perimental results are consistent with the expansion dynamics of an initially 
compressed nucleus and a break-up density of 0.3po- 

Figs. [2] and || show another important difference between peripheral and 
central collisions. For high excitation energies both curves follow the T = 
^Ek in /A free gas limit shifted by the remaining binding energy. However at 
low energies the peripheral setup shows a phase-transition like behavior: with 
decreasing excitation energy the temperature stays constant, and at even lower 
energies the nuclear liquid relation is reproduced. 

We conclude that to a large extent the temperatures observed in multifrag- 
mentation can be understood from the expansion dynamics into the spinodal 
regime and across the freeze-out density both for peripheral and for central 
collisions, despite the different break-up mechanism. The observed "phase- 
transition" plateau depends on the initial compression of the remnant, and for 
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usual central collision is missing. 
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